Bulk motion of granular matter in an agitated cylindrical bed 
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Experimental results are reported for the bulk motion induced in a bed of granular matter con- 
tained in a cylindrical pan with a flat bottom subjected to simultaneous vertical and horizontal 
vibrations. The motion in space of the moving pan is quantified. A number of distinct bulk dynami- 
cal modes are observed in which the particle bed adopts different shapes and motions. At the lowest 
pan excitation frequency ui, the bed forms a "heap," and rotates about the cylinder axis. As ui is 
increased, a more complex "toroidal" mode appears in which the bed takes the shape of a torus; in 
this mode, circulation occurs both about the cylinder axis, and also radially, with particles moving 
from the outer edge of the pan to the centre on the top surface of the bed, and back to the outer 
edge along the pan bottom. At the highest to, surface modulations ("surface waves" and "sectors") 
of the toroidal mode occur. The origin of this family of behavior in terms of the pan motion is 
discussed. 
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I. INTRODUCTION 

Granular materials play an important role in many in- 
dustrial processes. In applications, granular materials 
are often subjected to vibrations to generate material 
transport, mixing, or size segregation. Consequently, the 
fundamental properties and motions of granular matter 
confined to a vibrating container have been widely stud- 
ied. Many fundamental studies have focused on the mo- 
tion of particle beds in a rectangular container (the pan). 
Typically, the bed is excited by pan vibrations along 
a single axis. Reviews of such studies can be found in 
Refs. [Q, ||, ||, [|, ||, 0]. Agitated beds of granular matter 
have been observed to adopt a wide spectrum of shapes 
and motions, such as heaping, convection, small ampli- 
tude wave motion, arching, and large amplitude wave 
motion; see e.g. Ref. R. These modes certainly depend 
on the nature of the vibrations. However, their depen- 
dence on the pan geometry is less well understood. For 
example, in the case of purely vertical vibrations, the ef- 
fect of the pan shape is absent or weak [||, but this is 
unlikely to be generally true in the presence of horizontal 
vibrations. 

The motions observed in a cylindrical pan subjected 
to simultaneous horizontal and vertical vibrations, are 
reported in the present study. There are two principal 
motivations for studying this case. First, the response of 
a bed of granular matter to excitations occurring simul- 
taneously along three Cartesian axes has not been ex- 
tensively studied and remains poorly understood. How- 
ever, in many industrial applications, tri-axial excita- 
tions are highly desirable to optimize process parame- 
ters. Many fundamental studies have been made of sys- 
tems subjected to purely vertical or horizontal vibrations, 
but in only a few cases are vibrations in both directions 
addressed. An example is the work of Tennekoon and 
Behringer |^], who reported that when granular matter 
is subjected to simultaneous horizontal and vertical sinu- 



soidal vibrations, the phase difference between the com- 
ponents of vibration in the two directions becomes a key 
control parameter for the resulting motion. A quantita- 
tive understanding of the effect of such combined vibra- 
tions is therefore important for predicting and controlling 
the behavior of agitated beds in real applications. 

A second motivation is to study the motion generated 
in a container of cylindrical symmetry, with the cylinder 
axis oriented vertically. Again, this geometry is common 
in industrial devices such as sifting machinery, yet is less 
well studied in a research context compared to square 
and rectangular geometries. Note that the present case 
is distinct, both in motion and in orientation, from the 
more widely studied case of a continuously rotating drum 
having the cylinder axis oriented horizontally. 

A cylindrical geometry provides the opportunity to 
study phenomena unlikely to occur in the rectangular 
case. In particular, the absence of discontinuities (i.e. 
sharp corners) on the vertical boundaries of a cylindri- 
cal pan makes possible smooth circulatory motions of the 
agitated bed. These circulation modes are of interest for 
applications involving sifting and mixing, and for devel- 
oping agitated bed devices for continuous, as opposed 
to batch, processing. Initial studies have also demon- 
strated that these modes may display complex behavior. 
For example, Scherer et al. [Hj] have reported that when 
spheres are placed in a cylindrical pan and are subjected 
to horizontal shaking, they circulate in one direction at 
low packing densities, and in the other at high packing 
densities or high excitation frequencies. 



II. EXPERIMENTAL SET-UP 

The agitated bed apparatus studied here (Figs. [I] and 
|) is a modified SWECO finishing mill Model ZS30S66, 
a commercially available device widely used in industrial 
applications for sieving and polishing. It consists of a 
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(a) 



FIG. 1: Photographs of the agitated bed apparatus, viewed 
from above: (a) without particles; (b) with particles. 



rigid, circular, stainless steel pan of radius R = 0.381 m 
with a flat bottom. The pan "floats" on an array of nine 
springs spaced equidistantly around the circumference of 
the pan bottom. Vibrations are excited by a 0.5 HP 
motor mounted below the center of the pan. 

In order to induce simultaneous horizontal and verti- 
cal vibrations, and concurrently obtain control over their 
relative magnitudes, metal weights are mounted eccen- 
trically at the top and bottom ends of a shaft running 
through the center of the motor. The motor housing is 
rigidly attached to the pan via a flange, while the shaft 
rotates freely with respect to the pan. When the shaft ro- 
tates, the effect of the weights is to impose an oscillatory 
torque that tends to make the motor shaft deviate from 
the vertical. The affect of this torque is to cause the 
center of the pan to revolve about a vertical axis, and 
simultaneously cause the bottom of the pan to deviate 
from the horizontal. A complete description of the pan 
motion is given in the next section. The angle between 
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FIG. 2: (a) A schematic view of the apparatus. (b) 
Schematic representation of the relative positions of the motor 
and top and bottom weights, and the definition of (j>. 



the top and bottom weights, referred to here as the "lead 
angle" 0, can be varied from 0° to 180° [Fig. |(b)]. The 
motor shaft rotates the weights in a range of frequency 
lo of 10 to 20 Hz. 

The bottom surface of the pan is a stiff, stainless steel 
seive. The particles used in this study are too large to 
fall through the holes of the seive. The springs that bear 
the load of the entire pan-motor assembly are mounted 
on a rigid basement. Note also that there is a cylindrical 
"hub" attached to the center of the pan bottom, of radius 

10.4 cm and height 3.0 cm. 

The masses of the top and bottom weights attached to 
the motor are fixed in this study. As will be shown in 
the next section, the angle <fi between the weights controls 
the ratio of the amplitudes of the horizontal and vertical 
vibrations imparted to the pan. The rotational frequency 
lu of the motor controls the amount of energy introduced 
into the system. The principal control parameters for the 
experiments described here are thus <p and oj. 

The experiments are conducted with two sets parti- 
cles: "pill-shaped" oblate spheroid particles, and spheri- 
cal particles. The oblate spheroids have a major axis of 

13.5 mm and a minor axis of 7 mm, a mass of 0.9 g and a 
density of 1200 kg/m 3 , with a hard and smooth polished 
surface. The spheres have the same volume as the oblate 
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spheroids, a density of 1100 kg/m 3 , and a hard, smooth 
but unpolished surface. Although the qualitative behav- 
ior observed in both cases is similar, it is found that the 
shape of the particles influences some quantitative prop- 
erties of the motion. Thus the results obtained with each 
type of particle are presented separately. 

The mass of the unloaded pan is approximately 150 kg. 
All tests are conducted for particle bed loads below 45 kg 
and it is found that the influence of the bed load on the 
pan motion is negligible. Also, it is confirmed in all tests 
that the pan moves as a rigid body. 

The motion of the pan is monitored using Entran ac- 
celerometers: one tri-axial (model EGA3-F-10-/5) and 
three mono-axial (model EGA-F-10-/5) units are used. 
The location of the tri-axial accelerometer is shown in 
Fig. |^; the others are positioned at various locations 
along the top edge of the pan, according to the degree 
of freedom to be measured. The accelerometer data are 
collected using a computer-based data acquisition sys- 
tem and processed to obtain the velocity, displacement 
and phase-shift data for the pan motion. A sampling 
frequency of 2 kHz is used for data acquisition, allow- 
ing a minimum resolution of 165 points/cycle. This re- 
sults in an error in the phase angle measurements of ±3°. 
The uncertainty in the acceleration and bed displacement 
are estimated to be ±0.1 g and ±0.05 mm, respectively, 
where g is the acceleration due to gravity. 

The surface motion of the particle bed is studied by 
visually tracking particles. Differently colored particles 
are used to obtain information on the average mixing 
time. For the most interesting modes of motion, a lo- 
calized group of particles is coated with a fluorescent dye 
and the subsequent motion of this group is recorded pho- 
tographically. The dye is sodium fiuorescene, which is 
excited using an ultra-violet light source and observed in 
the 532 nm visual range. 

The results of the experiments are presented in two 
stages. First, a quantification of the pan motion is pre- 
sented in Section 3. Then the modes of particle motion 
that are observed in the bed are described in Section 4. 



III. PAN MOTION 

As described above, the experimental apparatus is a 
commercially available unit that generates simultaneous 
horizontal and vertical vibrations of the pan. The first 
task of the present study is therefore to quantify the mo- 
tion generated by this device. To achieve this, two coor- 
dinate systems are defined: one fixed with respect to the 
pan (the "pan frame" ) , and one fixed with respect to the 
laboratory (the "lab frame"). 

The intersection of the cylinder axis of the pan with 
the plane defined the by top edge of the pan is chosen as 
the origin of the Cartesian coordinate system [x p , y p , z p ) 
fixed in the pan frame. The cylindrical symmetry of the 
pan is broken by the existence of a spout located at a 
point on the circumference of the pan. The x p -axis in the 
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FIG. 3: Photographs of laser traces visualizing the circular 
motion of the pan center for (a) <f> — 10°, U) = 30 Hz; (b) 
(j> = 60°, lo = 30 Hz; and (c) cj> = 100°, u = 30 Hz. 



pan frame is thus defined as the line that passes through 
the origin and the location on the spout where the tri- 
axial accelerometer is mounted. The y p -axis is the line 
perpendicular to the a; p -axis in the plane defined by the 
top edge of the pan, and the z p -axis is perpendicular to 
both the x p and y p -axes. The Cartesian coordinate sys- 
tem (x, y, z) in the lab frame is defined as that which is 
coincident with the pan frame when the apparatus is at 
rest. 

The accelerometers provide data that are analyzed to 
give the lab frame coordinates as a function of time for 
points fixed in the pan frame. To specify the motion of 
the pan, in the following attention is restricted to points 
fixed in the pan frame lying in the z p = plane, since 
this is the plane in which the accelerometers are located. 
Based on the accelerometer data, the motion in the lab 
frame of such a point is found be consistent, within the 
error of measurement, with the following model: 



x(t) 

y(t) 

z(t) 



where, 



and 



x p + x max cos(c^) 

y P + 2/max COs(ujt + a xy ) 

A(r p )cos(tut + a xz - 6), 



Mr P ) = £ 



(1) 
(2) 
(3) 



(4) 



(5) 



In the above equations, x max and y max are the ampli- 
tudes of the x and y displacements of the pan-frame ori- 
gin in the lab frame; a xy is the phase shift between the x 
and y displacements in the lab frame; and 9 is the angle 
measured about the origin from the a; p -axis to the point 
(x p ,y p ). a xz is the phase shift between the x and z dis- 
placements in the lab frame, as measured on the a; p -axis; 
thus (a xz — 9) is the phase shift between the x and z 
displacements in the lab frame, to be found at the point 
(x p ,y p ). A(r p ) is the amplitude of the z displacement of 
a point in the z p = plane that is a distance r p from the 
origin in the pan frame. R — 0.381 m is the distance from 
the origin in the pan frame to the edge of the pan, and 
z max is the amplitude of the z displacement of a point on 
the edge of the pan (i.e. at r p — R). 
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FIG. 4: Phase shift function of (f> for various w 



FIG. 6: Phase shift a xz , as a function of <j> for various ui. 
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FIG. 5: Maximum displacements a; ma x, 2/max and z max as a 
function of c/>. These data are averages over approximately 20 
runs over the range of ui studied. No systematic variation is 
observed as a function of uj, though individual measurements 
are scattered about the average by about 15%. 



The form of Eqs. [I] and is motivated by the obser- 
vation that under all operating conditions, the center of 
the pan appears to execute an approximately circular 
orbit in a horizontal plane. To visualize this, a laser 
source is mounted at the center of the pan, pointing up- 
ward, normal to the pan bottom. The motion of the pan 
causes the laser to trace out a pattern on a screen posi- 
tioned horizontally above the apparatus. Time-exposure 
photographs of these traces are shown in Fig. || for vari- 
ous operating conditions, and are consistent with circular 
motion of the pan center. 



If the motion is nearly circular, the magnitude of the 
phase shift a xy should be approximately 90°. This is 
confirmed in Fig. ^, where a xy is plotted as a function 
of <p for various u). An apparatus with perfect cylindrical 
symmetry would have a xy = 90°. Small asymmetries 
in the apparatus, the most important of which is likely 
the spout on the edge of the pan, result in the observed 
variation in the range 93° < a xy < 100°. 

Circular motion of the pan center also requires a; max 
and 2/ m ax to be equal. Measurements of £ max and i/ max 
at r p = R (Fig. ||) show that this is observed to be true 
within the measurement error across the range of <p stud- 
ied here. Note that as u> varies in different test runs, x max 
and y max are found to be scattered around the average 
values reported in Fig. ||by about 15%; however, no sys- 
tematic variation with u> can be identified. Tests also 
show that the values of x max and y m ax are independent 
of the position around the edge of the pan at r p = R at 
which they are measured, consistent with the description 
in Eqs. [j] and g. Note that for x max or j/ max to have the 
same measured value independent of r p at the point of 
measurement requires that the tilt angle of the z p — 
plane with respect to the z = plane be small. This is 
confirmed below. 

Eq. |^ characterizes the tilting motion of the pan by 
quantifying the vertical deviation of a point fixed in the 
pan frame (in the z p = plane) with respect to the z = 
plane in the lab frame. As in the case of the x and y mo- 
tions, the accelerometer data for z(t) indicate a sinusoidal 
function of t. The form of Eq. |^ also assumes that the ori- 
gin of the pan frame does not leave the z — plane of the 
lab frame; this assumption is confirmed by direct visual 
observation using a high-speed (500 frames/s) camera. 
From geometric considerations, the amplitude A of the 
z variation in Eq. |3] is a function of the distance r p of 
the point of measurement from the center of the pan, as 
indicated in Eq. 13. The phase shift between the x and 
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FIG. 7: Schematic representation of the pan motion. The 
dashed circle represents the orbit of the pan centre in the lab 
frame. The solid circle represents the rim of the pan. 



z displacements is a difference of two contributions, a xz 
and 9. a xz is the phase shift between the x and z dis- 
placements as measured on the x p axis at r p = R. 9 is 
a geometrical term that accounts for the difference be- 
tween the phase shift observed on the and that 
which would be observed at an arbitrary point [x p , y p ) 
in the z p = plane. As shown in Fig. |6|, a xz does not 
depend strongly on oj, but does depend on <f>. 

The net effect of the motions described by Eqs. (1-3) 
are summarized in Fig. [?]. The motion of the pan can be 
thought of as a superposition of two motions: 

(i) The center of the pan revolves about the lab frame 
origin in a nearly circular, horizontal orbit. This motion 
is directly described by Eqs. (1-2). 

(ii) The pan is tilted at an angle 7 in such a way that 
a unit normal vector rooted at the pan center precesses 
about the vertical at the same frequency at which the cen- 
ter orbits the origin, but phase shifted with respect to the 
orbital motion of the center. The fact that 7 maintains 
a constant value can be shown from Eq. (3) by consid- 
ering (e.g.) the case r p = R and setting 9 = ujt + a xz . 
These conditions describe the point on the edge of the 
pan that at any given t has the largest positive z displace- 
ment. For all t, this displacement is a constant, z max , 
and so 7 = sin _1 (z max /i?) is also constant. Note that for 
the apparatus and operating conditions described here, 
R = 0.381 m, and z max is never greater than 0.0035 m, 
giving 7 < 0.5°. This confirms the assumption of small 
tilt angle required for the chosen form of Eqs. (1-2), as 
indicated above. 

The measurements shown in Figs. show that both 
the amplitudes and phase shifts required to specify the 
pan motion using Eqs. (1-3) are mainly controlled by 
4>. The effect of u> is less problematic, since the ratio of 
displacements £ max and z max do not vary in a systematic 
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FIG. 8: Diagram identifying the regions of the u-<j> plane 
in which distinct bulk dynamical modes of the particle bed 
are observed for (a) oblate spheroids, and (b) spheres. Note 
that in (a) a region of dynamic equilibrium is observed near 
cj> = 10°, where the bed does not perform any net rotation 
about the pan centre. This is the boundary region between 
counter-clockwise and clockwise bed rotation ("laps"). 



way with ui. However, since w controls the rate at which 
energy is introduced into the system, both (f> and u> will 
play a major role in determining the particle motion, as 
shown next. 



IV. MODES OF PARTICLE BED MOTION 

A diagram showing the observed bulk dynamical 
modes of the particle bed as a function of uj and <f> 
is shown in Fig. |§|(a) for the oblate spheroids, and in 
Fig. |^(b) for spheres. The bulk motion of the bed is a 
strong function of the ui. Four distinct modes of bed 
motion are observed, here termed "heaping," "toroidal 



FIG. 9: Fluorescene dye visualization of the particle bed sur- 
face motion in the heaping mode. 30 kg, <j> — 60°, uj — 30 Hz. 



motion," "surface waves" and "sectors." Each of these is 
described in detail below. The range of the uj at which 
each of the modes appears depends on the particle shape 
as can be seen from Fig. but the qualitative behavior 
of the modes is mostly similar. In the following sections, 
the motion in general is discussed, and differences due to 
particle shape are identified where appropriate. 



A. Heaping 

At the lowest uj studied the particles form a heap with 
an elevated center. The relative particle movement and 
mixing is very small and the particle bed rotates almost 
as a solid body. This motion is visualized in Fig. ^|. Flu- 
orescent dye is continuously added at a specific location 
to the surface of the particle bed while in steady-state 
motion. The subsequent motion of the dyed particles is 
then captured photographically under ultra-violet light, 
as shown in Fig. |j| The dyed particles on average un- 
dergo a slow processional motion about the pan center. 
The radial component of the particle motion and the dif- 
fusion of particles is negligible. Notably, for most 10 and 
4>, the particle bed rotates in the direction opposite to 
that of the orbital motion of the center of the pan; the 
exceptions to this occur at the lowest 4> an d uj studied 
(e.g. cf> — 10° and 16 Hz for oblate spheroids), where the 
rotational directions are the same. The direction of mo- 
tion of the bed appears to depend on the phase difference 
a yz = olxz — ct xy . For a yz > 0°, the bed rotates in the 
direction opposite of the pan center. 

The slope of the particle bed (from the center to the 
outer edge) increases with an increase in <j>. For example, 



FIG. 10: Fluorescene dye visualization of the particle bed 
surface motion in the toroidal mode. Particles spiral to the 
center of the bed along the surface and toward the edges along 
the pan surface. These reappear at the surface of the bed 
along the edges. 30 kg, <f> = 60°, uj = 56 Hz. 

the bed of oblate spheroids attains the maximum slope 
for <j) = 100°, and a nearly flat surface when cf> = 10°. 
Fig. H schematically shows the cross-sectional shape of 
the surface during the heaping mode. 

The spherical particles move in more tightly packed 
layers than the oblate spheroids. For = 100°, the layers 
are very tightly packed and the local arrangement of par- 
ticles is similar to a close-packed lattice. As 4> increases, 
the uniformity of the lattice breaks down. The bed is 
still tightly packed at <fi = 100°, but more pronounced 
relative motion of the particles is observed. This change 
in behavior may be related to the increase of z m ax as <j) 
increases (see Fig. ||). 

Previous studies of granular matter in vertically ex- 
cited rectangular containers |l2| have shown that rela- 
tive particle motion is observed when the maximum ver- 
tical acceleration exceeds 1.2g. In the present tests, it is 
found that the local acceleration is below this threshold 
for uj < 10 Hz. However, at <f> = 100°, the acceleration 
at the edge of the pan does approach 1.2g. 



B. Toroidal Motion 

As uj increases the particle bed undergoes a transfor- 
mation to a toroidal shape and the particle motion be- 
comes organized. This motion is highly coherent as seen 
in the fluorescent dye visualizations in Fig. |l^. Dye is 
continuously injected onto the surface of the moving bed. 
Particles initially on the outer edge move in a spiral mo- 
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FIG. 11: Schematic of cross-section of particle bed, showing 
direction of particle flows in toroidal motion. 



tion to the center. Near the center they are subducted 
and travel outward along the bottom of the pan, to be 
re-entrained along the walls and reappear at the surface 
at the outer edge of the pan. As seen for the heap, the 
bulk rotational motion of the bed in the horizontal plane 
is opposite to that of the center of the pan. The entire 
particle bed adopts the shape of a torus, through which 
individual particles move in helical trajectories. The par- 
ticle motion in a cross-sectional plane through the toroid 
is depicted schematically in Fig. |ll|. Note that in addition 
to the primary circulation roll occupying the center por- 
tion of the bed, a small secondary (and counter-rotating) 
circulation roll forms on the top of the primary one, near 
the pan wall. The particles move toward the center along 
the top surface in the primary roll, and toward the wall 
along the top surface in the secondary roll. The slope 
of the particle bed surface is greater in magnitude, and 
opposite in sign, than that observed for the heap motion. 

In this mode of motion, the particle bed is fluidized. 
The increased agitation results in a relatively rapid dif- 
fusion of the particles and a net increase in the mixing 
rate. 

The toroidal motion observed for the circular pan has 
not been observed in rectangular geometries. An inward 
spiraling motion has been observed for granular media 
underneath rotating fluids Jl(| and the radial segrega- 
tion of granular mixtures has been observed in rotating 
cylinders For vertically excited rectangular beds, 

stationary convection cells, similar to Rayleigh-Benard 
instabilities, are observed |l3). For horizontal vibra- 
tions, experimental 14 and computational fHI studies 



have shown that the convective motion consists of gran- 
ular matter rising to the surface at the middle of the 
container. In contrast, in the present study, the particles 
rise along the walls and are subducted at the middle of 
the pan. 




FIG. 12: Two fluorescene dye visualizations of sectors. The 
sectors can be visualized because the particles located on the 
crests dry faster than those in the trough and loose fluores- 
cence. 



C. Surface Waves 

For the two modes described above, the bed height is 
only a function of the distance from the center of the 
pan, and the shape of the bed is symmetric about the 
cylinder axis of the pan. However, for uj > 15 Hz, crests 
and troughs appear on the surface of the bed as a function 
angular direction around the pan. The crests, or fronts, 
propagate in the same direction but at a higher speed 



than the bulk rotation of the bed in the horizontal plane. 
The speed of the particle fronts initially increases as uj 
increases, but then decreases as uj approaches 18 Hz. The 
amplitude (i.e. the height of the crest above the bed) 
decreases as uj increases. Unfortunately, this motion is 
very difficult to capture on still media. 

The individual particle trajectories can still be de- 
scribed as helical. The motion of the crests appears to 
be similar to a traveling wave. 



8 



D. Sectors 

For 18 Hz< u> < 20 Hz, the motion of crests on the 
surface of the bed becomes nearly stationary, dividing 
the bed into well-defined "sectors." These sectors are 
visualized using fluorescent dye in Fig. [l^. The number 
of sectors present in the bed is found to be a function of 
4> and u>. The sectors are stable and easily recognized for 
the oblate spheroids. However, for the spherical particles, 
stationary sectors are more difficult to achieve and could 
not be reliably reproduced. 

Although the general bulk motion of the particles 
remains helical, the sector mode is distinct from the 
toroidal mode. Defining the general direction of the pro- 
cessional motion as streamwise, it is seen that the parti- 
cles on the upstream side of the crests are drawn into the 
bed, while on the downstream side particles emerge on 
the surface. The upstream particles have been exposed 
to air for a longer period and have dried, causing the 
dye to loose florescence. The upstream side of the crests 
thus appear darker than the downstream side, yielding 
the "rays" seen in Fig. |l2|. Also, in the sector mode, the 
secondary circulation roll near the wall that was charac- 
teristic of the unmodified toroidal mode of motion could 
not be observed. 

As uj is further increased to 26 Hz, the toroidal motion 
(without sectors or surface waves) is recovered. Unfortu- 
nately, the mechanical limitations of the apparatus did 
not allow testing at larger lo. 

V. DISCUSSION 

A. Rotation of the particle bed in the horizontal 
plane 

The direction of rotation of the particle bed in the 
horizontal plane is generally opposite to that of the center 
of the pan. This behavior is initially counter-intuitive, 
and explaining it requires a detailed understanding of the 
interplay between the bed motion, the pan orientation, 
and the horizontal and vertical components of the pan 
motion. 

A crucial element for resolving this behavior is to de- 
termine the shape and motion of the zone of contact of 
the particle bed with the pan. Because the pan is tilted 
away from the vertical at all times during its motion, it 
seems likely that only a particular sector of the bed is in 
contact with the pan at any given time; and, that the an- 
gle at which maximum contact occurs between the bed 
and the pan rotates about the pan center at the same 
frequency oj as the pan center itself. If this is true, at 
least two contributions to the mechanism of bed rotation 
can be envisioned: 

(i) Deflection by the pan surface: Consider an element 
of the pan surface within the zone of contact between the 
pan and bed. At any given time, the normal to this suface 
element can be decomposed into vertical, radial, and tan- 




FIG. 13: Schematic illustration of scenarios for explaining 
direction of bed rotation. The large circle represents the edge 
of the pan as viewed from above, c is the circular trajectory 
of the centre of the pan as seen in the lab frame, while d 
is the instantaneous direction of motion of the pan centre 
at the moment depicted, s is a surface element of the pan 
bottom along the direction z of the zone of maximum contact 
between the pan and the bed. Scenario (i): If the vector t 
represents the tangential component of the normal to s, then 
particles will be deflected along t, leading to bed rotation in 
the direction b, opposite to c. Scenario (ii): If the vector t 
represents the tangential component of the velocity of s, then 
particles will be "dragged" in the direction of t, leading to 
bed rotation along b, again opposite to c. 



gential components, defined with respect to a cylindrical 
coordinate system in the lab frame. Particles striking 
this surface element will tend to be deflected according 
the to the orientation of the surface normal. In the ab- 
sence of other mechanisms, a net rotational motion about 
the pan centre will occur if the direction of the tangen- 
tial component remains constant. This will be true if (as 
assumed above) the zone of particle-bed contact rotates 
around the pan bottom at the same frequency as the pan 
itself precesses. The direction of this bed rotation will be 
independent of the rotational direction of the pan center, 
since it depends on the (presently unknown) location of 
the region of particle-bed contact. Hence, the bed could 
well be set into a rotational motion opposite to that of 
the pan center (Fig. |ili| ). 

(ii) Entrainment by the pan surface: Independent of 
the above mechanism, the effect of friction between the 
particles and the pan bottom should also be considered. 
As described in Section ||, the pan bottom consists of a 
metal screen, and so is "rough" on the scale of the par- 
ticles themselves. Consider the same surface element in 
the contact zone as discussed above, but instead of its 
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FIG. 15: A plot of the overs to lap ratio for the oblate 
spheroids. 



orientation, consider its velocity vector decomposed into 
vertical, radial, and tangential components. This rough 
surface clement will tend to transfer momemtum to par- 
ticles in contact with it according to the direction of its 
velocity. Analogous to the reasoning given above, parti- 
cles in the contact zone will be subjected to a constant 
tangential surface velocity, yielding a net rotation of the 
bed about the pan center. Also as above, the direction of 
the resulting bed rotation will be independent of that of 
the pan center, and so it should be possible to establish 
conditions where these rotations would be in opposite in 
direction (Fig. |l3|). 

Whatever the underlying mechanism, an explanation 
can be tested against the observed dependence of the mo- 
tion on the control parameters. For example, the speed 
of rotation of the bed in the horizontal plane depends on 
the phase difference between the horizontal and vertical 
pan motions, a xz . For 0° < a xz < 180°, the particles 
are observed to move in the direction opposite to that of 
the bed. However, as mentioned earlier, it is indeed pos- 
sible to realize conditions where the bed rotates in the 
same direction as the pan, when <f> is changed such that 
a xz < 0°. 

In addition, the rotational speed of the bed in the hor- 
izontal plane increases with lj. Furthermore, as shown in 
Fig. |l4|, the rotational speed is found to be greatest when 
ctxz = 90°, corresponding to = 60° (see Fig. |J). 

Elucidation of these phenomena requires information 
on the (time-dependent) location of contact between the 
pan and the bed. This information is not available in 
the present measurements. These questions are being 
explored through further experiments, and through com- 
puter modelling. 



B. Radial motion of the particles 

Understanding the radial motion of the particle bed 
also requires a careful consideration of the bed and pan 
motions. In the toroidal mode, the particles move along 
the bottom of the bed toward the wall, and then back 
toward the center of the pan along the top surface of the 
bed. It seems likely that the motion of particles near 
the pan bottom is due to their being entrained by the 
motion of the pan surface itself. The inward motion of 
particles on the bed surface, where the particle packing 
is looser, may be dominated by gravity-induced down- 
ward flow toward the bed center, where the bed depth 
is smallest. However, these possible explanations do not 
fully elucidate the origin of the toroidal shape adopted 
by the bed in this mode, and so further study of this 
behavior is needed. 

The ratio of the rotational speed of the particle bed 
in the horizontal plane, to the radial speed is show in 
Fig. |l5|. To measure this, the motion of color-tagged 
particles is observed in steady state. The average time 
tl required for a tagged particle to complete one circuit 
(a "lap" ) around the pan in the horizontal plane is mea- 
sured; the average time to required for a tagged particle 
to complete one circuit (an "over") from the center of 
the pan to the outside edge and back to the center is also 
measured. The ratio tl/to is shown in Fig. [l5| 

In the heaping mode the radial motion is negligible, 
as indicated in Fig. ^. When the vertical acceleration is 
sufficiently high to fluidize the particles (i.e. the particles 
are not in continuous contact with the pan surface), the 
radial and horizontal rotational speeds of the particles 
are found to depend only on the phase difference between 
the horizontal and vertical pan motions. This is shown in 
Fig. [L5f as bj is increased, tq/tl approaches a constant 
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g as a function of if). 

value, which is a function of <fi. 

C. Critical acceleration and fluidization 

The data in Fig. Hare used to determine the critical 
radius, shown in Fig. |16|, at which the maximum vertical 
acceleration exceeds g. For <fi = 10° and u> < 18 Hz, the 
vertical acceleration is less than g and the particles at 
the bottom of the bed are constantly in contact with the 
pan surface. This is consistent with the radial motion of 
the particle bed being negligible, as indicated in Fig. 0. 
Based on the model motion in Eqs. |l[|3[ for lo = 20 Hz, 
the vertical acceleration exceeds g when r p /R > 0.81. 
Under these conditions the toroidal motion is induced, in 
which case an overs-to-lap ratio could be clearly defined 
(Fig. |l5]). As </> increases, the critical value of g occurs in 
the bed at smaller u>, in agreement with the observations 
summarized in Fig. |§|. 

When the maximum vertical acceleration in the pan is 
less than g, the tilting motion of the pan induces only a 



rotational motion of the bed in the horizontal plane, since 
the particles at the bed bottom remain more or less in 
constant contact with the pan surface and the net radial 
bed motion is zero. When the maximum acceleration ex- 
ceeds g, the particles at the bottom of the bed are lifted 
above the pan surface, allowing a net radial displace- 
ment to be induced, and leading to the toroidal motion. 
Note that the toroidal mode is induced below the critical 
acceleration of 1.2g reported for fluidization of granular 
matter in other studies f6|, 0, [if]. For <j> > 30° and 
u) > 16 Hz, the particle bed is fluidized over most of the 
pan area. This is the range of uj in which the "wave" and 
"sector" modes are seen. 

In their study of a deep bed of granular matter, Wass- 
gren et al. have shown that the bed undergoes a stage 
with an arched base, and during this stage, the bed moves 
out-of-phase with the container. A similar argument can 
be put forth for the circular bed in the present study. 
The motion that is imparted to the bed depends on the 
pan motion and the phase difference between the pan 
and the bed. When the bed moves out-of-phase with the 
pan, the point of next contact with the pan becomes a 
function of the pan motion. The instantaneous traces of 
the pan vertical displacement show period-doubling when 
the surface waves were observed in the particle bed. At 
every second cycle, the maximum upward vertical dis- 
placement is slightly less than that of the first cycle. The 
difference in amplitude decreases as the surface waves ac- 
celerate and the sectors eventually appear. The sectors 
may thus correspond to a standing wave. 

Our observations are thus consistent with the existence 
of well defined states based on the frequency of the ex- 
citation. Computer simulations are currently underway 
to further elucidate the origin of the behavior described 
here 
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